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INTRODUCTION 

With  on*  notabl*  exception,  this  year's  activities  on  Grant 
APOSR-80-0282  were  much  different  than  the  activities  of  preceeding 
years.  This  fact  is  a  simple  consequence  of  the  evolutionary  nature  of 
this  program.  The  first  few  years  were  spent  demonstrating  the 
feasibility  of  the  experiment  and  recording  and  evaluating  velocity 
and  pressure  data  to  be  used  as  the  forcing  function  for  the  boundary 
layer  itself 1 »3.  During  the  next  year  the  initial  boundary  layer  data 
were  recorded  for  a  variety  of  test  conditions3.  This  was  a  very 
useful  and  exciting  period,  since  much  of  the  data  represented  a  new 
combination  of  parameters.  The  last  year  was  spent  largely  in 
consolidating  the  data,  "filling  in  the  envelope,”  as  it  were.  It  was 
also  a  period  of  reflection,  consolidation  and  improved  understanding  of 
both  the  older  and  the  newer  data,  which  included  studies  of  additional 
features  of  the  experiment,  such  as  near  wake  velocity  profiles  and 
separated  flow  pressure  distributions.  Additional  tests  were  also 
conducted  to  help  delineate  the  velocity  part  of  the  unsteady  forcing 
function  from  the  acoustic  or  pressure  part.  In  some  circumstances 
tests  were  conducted  either  to  demonstrate,  say,  the  two-dimensionality 
or  lack  of  it,  on  the  metric  surface  or  to  confirm  or  reject  hypotheses 
or  explanations  of  various  features  i,n  the  data.  These  results  will  be 
discussed  briefly  below.  Some  interpretations  of  this  data  will  also  be 
briefly  set  forth. 

The  notable  exception  referred  to  in  the  opening  sentence  is  a 
result  of  a  failure  to  either  find  a  suitable  research  student  or  to 


t*  '•  i  i. . ;  - 1  , 

"hief,  Te 


bMU'al  Information  Divisi 


hire  a  suitable  engineer  to  complete  the  development  of  the  floating 


element  shear  gauge.  The  experience  from  the  proceeding  two  years 
suggested  that  unless  a  suitable  student  or  engineer  was  available,  the 
results  were  not  worth  the  cost.  This  short  term  loss  now  appears  to 
have  resulted  in  a  long  term  gain  for  two  reasons.  First  a  cheaper, 
more  reliable  method  has  been  adapted  for  our  use,  and  second  we  have 
reprogrammed  the  assets  and  improved  our  ability  to  take  simultaneous 
data  from  many  hot  wire  and  pressure  sensors.  Such  data  appear  to  be 
very  valuable  in  the  region  of  near  separation  in  unsteady  flow. 

In  this  report  the  important  aspects  of  the  last  year's  studies 
will  be  briefly  outlined.  The  results  of  the  last  few  years  will  be 
placed  in  a  more  general  context.  This  discussion  will  lead  to  a 
discussion  of  some  important  unanswered  questions  and  their  importance. 
Technical  Summary 

The  important  technical  points  to  be  discussed  are  as  follows:* 

1.  The  nature  of  the  various  correlation  terms  in  unsteady  flow. 

2.  Experimental  demonstrations  of  the  two  dimensionality  of  the  flow 
field. 

3.  Experimental  delineation  of  acoustic  effects  from  unsteady  velocity 
effects. 

4.  Characteristics  of  unsteady  separated  flow  in  this  experiment. 

♦The  geometry  and  the  detailed  nature  of  the  experimental  conditions 
and  data  acquisition  system  are  described  in  detail  in  References  2 
and  3. 
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5.  Consents  on  the  trailing  edge  flow. 

6.  On  the  proposed  method  of  measuring  shear. 

Results  of  a  study  of  the  turbulent  unsteady  near  wake  completed  in  the 
fall  of  1982  are  given  in  the  attached  copy  of  AIAA  paper  83-0127. 


On  the  nature  of  the  several  correlation  terms 

Recall  that  our  signals  for,  say,  the  tangential  component  of 
velocity  are  defined  to  contain  three  parts  (Refs.  4,  5,  6):  a  time 
mean  part,  a  periodic  but  not  necessarily  sinusoidal  part,  and  the  part 
left  over.  Thus  if  u(x,t)  is  our  measured  signal,  then  the  time  mean  is 
defined 

i  rT 

u(x  )  ■  lim —  I  u(x,t)dt 

T-*-co  21  J  _<p 

Further  we  can  take  an  ensemble  or  phase  lock  average 

1  N_1  2 

<u(x,t)>  ■  —  2  u  *x»t  +  n  ~  +  9) 

M  n-0  u 

where  0  9  <_  2w/u  ,  but  is  constant  for  each  average.  Note  then 

<u(x,t)>  -  u(x)  +  u ( x , t ) 

This  equation  is,  in  fact,  the  definition  of  the  periodic  part  of 
the  signal,  u.  The  remaining  term  is,  by  definition 

u'(x,t)  “  u  -  <u>  -  u  -  u  -  u 


Note  u1  has  two  properties,  by  definition 


A) 


u'(x,t)dt  -  0 


B) 


i  rT 

lim  —  i 

T+m  2T  J -T 

N-1 

Y  U'  (x,t  +  n  —  +  0) 
n-0  “ 


That  la  to  say  u'  has  zaro  time  mean  and  has  zero  ensemble  average. 

We  are  assuming  that  u1  corresponds  to  the  usual  turbulent  velocity 
component  that  one  measures  in  steady  flow.  In  turbulent  boundary 
layer  theory  the  term  <uv>  is  of  great  importance  since  it  leads  to  the 
"Reynolds  Stress”  terms.  Based  upon  the  decomposition  above  we  find 


IS 


<uv>  ■  uv  +  uv  +  uv  +  uv  +<uv  >  +  <uv  >  +  <uv  > 
+  <uv  >  +  <u  V  > 


s 


i 


The  time  average  of  this  term  gives 


lim  — 
m _  2T 


/T  _ *  i  f  T 

(uv)dt  -  uv  +  lim  —  I  (uv  +  uv  +  uv)dt 

-T  T-*o»  2T  J 


+  lim 
T-*oo 


—  r 

21  J -] 


(uv  +  u  v  +  u  v  )dt 


•tc. ).  The  terms  u'v*  is  the  classic  Reynolds  stress  term.  The 
last  two  terms  would  seem  to  have  a  small  value  (random  rectifications). 
The  real  issue  is  whether  or  not  u'  as  a  function  of  time  is  of  finite 
or  zero  measure.  In  the  latter  case  u'v  is  identically  zero.  In  the 
former  it  should  be  small.  According  to  our  data  terms  like  u'v  and 
uv'  are  less  than  10-4  0f  jjg  or  leS8  than  U«V' . 

Similarly  we  may  take  an  ensemble  average ,  subtract  out  the  mean 
and  end  up  with 

UV  *  UV  +  uv  +  UV  +  UV  +  U  V  +  U  V 
-  uv  ♦  uv  +  UV  +  U  V 

The  appearance  of  the  non-zero  term  u'v'  is  a  curiosity.  The  basic 
premise  is  that  the  periodic  part  of  each  term  in  u  and  v  had  been 
removed  by  the  previous  manipulation.  Figure  1  shows  typical  values  of 
across  the  wake  layer7.  It  is  definitely  not  zero.  It  was 
suggested  by  Lorber^  that  a  simple  geometric  argument  accounts  for 
this  ten.  To  follow  this  argument  keep  in  mind  that  u  is  really  a 
function  of  x  and  y  as  well  as  time.  In  the  unsteady  boundary  layer  we 
expect  each  of  the  several  thicknesses  characterizing  the  boundary  layer 
to  vary  over  a  cycle.  It  is  well  known  that  the  turbulent  intensity  is 
a  function  of  the  distance  from  the  wall.  Hence  the  probe  sees  a 
variation  in  scale  position  (l.e.,  y/6  with  y  fixed).  So  since  the  u' 
and  v'  correspond  to  different  heights  in  the  boundary  layer  over  a 
cycle,  so  u'v'  is  non-zero.  Thus,  there  is  an  unsteady  contribution  to 
the  Reynolds  stress,  at  a  fixed  y  that  varies  over  the  cycle.  This 


model  will  be  discussed  further  below. 

Experimental  Demonstrations  of  Two-Dlaenslonality 

The  two-dimensional  characteristics  of  the  flow  pattern  were 
Measured  in  two  ways.  First,  pressure  taps  were  mounted  in  the  lateral 

positions  ±12.5%  and  ±25%  of  chord  from  the  mid  span  location  at  6 

8 

chordwise  locations  between  x/c-0.025  and  0.98  .  Figures  2A,  B,  and  C 
show  the  comparison  of  this  data  for  attached  flow.  The  data  show  the 
flow  is  locally  two-dimensional.  Figures  3A-C  show  the  same  kind  of 
data  when  a  laminar  separation  bubble  is  present  near  the  leading  edge. 
This  separation  point  location  can  be  predicted  by  Stratford's 
criteria0  when  based  upon  mean  flow  properties  to  within  about  20%  (or 
about  1/2%  of  the  chordwise  position).  The  reattachment  zone  is 
predicted  by  Owen  and  Klanfer's  criteria10  equally  well.  Figures  4A-C 
show  results  for  massive  separation  of  an  artificially  tripped  boundary 
layer,  while  Figure  4D  shows  the  mean  pressure  distribution  with  natural 
transition.  Lack  of  two-dimensionality  in  the  leading  edge  bubble 
separation  has  been  observed  by  Winklaman  at  the  U.  of  Maryland11. 

More  important,  our  data  shows  the  separation  zone  for  the  case  of 
artificial  transition  is  more  three-dimensional  than  either  the 
laminar  bubble  or  the  massive  separated  region. 

The  other  procedure  to  test  the  two-dimensionality  is 
through  wake  measurements.  Typical  data  is  shown  in  Figures  5  and  6. 

The  conclusion  is  that  locally,  the  flow  is  two-dimensional  for  at  least 
the  distance  of  ±  one  quarter  span  on  either  side  of  the  airfoil 


centerline.  The  two  regions  selected,  the  location  of  the  separation 
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bubble  and  the  near  wake,  are  regions  where  one  would  expect  the 
departure  from  two-dimensional  flow  to  be  the  largest.  Hence  we  are 
confident  our  conclusion  that  the  flow  is  essentially  two-dimensional  is 
valid. 

Characteristics  of  Unsteady  Separated  Flow 

In  this  discussion  we  will  deal  with  bulk  characteristics  associated 
with  separated  flow.  This  restriction  is  primarily  due  to  problems 
associated  with  measuring  the  details  associated  with  unsteady  separated 
flow.  (It  is  likely  that  a  computer  contolled  laser  anemometer  could  track 
the  separation  point;  while  there  is  no  substantiation  of  this  conjecture 
on  separated  flaw,  such  measurements  have  been  made  in  the  wake  of  a 
propeller. )  Thus  separation  is  located  more  or  less  by  examination  of  the 
mean  and  fluctuating  pressure  distribution.  Hence  the  final  separation 
point  location  is  known  no  more  accurately  than  the  distance  between 
pressure  taps.  Figure  7  shows  the  variation  of  Cp  vs  x/c  as  a  function 
of  reduced  frequency,  k,  in  the  fully  attached  case8.  Note  the  high 
value  of  the  peak  suction  pressure  coefficient  (-6.2),  and  the  tendency  of 
the  mean  pressure  coefficient  difference  to  approach  zero  at  the  trailing 
edge.  Figure  8  shows  the  change  in  mean  difference  when  separation  is 
present  on  the  upper  surface8.  In  this  case,  the  mean  lower  surface 
pressure  coefficient  is  as  it  was  before,  except  near  the  leading  edge,  but 
the  mean  upper  surface  pressure  coefficient  is  radically  transformed. 
Figures  9  and  10  summarize  the  current  data  on  attached  flow,  separated 
flow  and  the  boundary  regions  dividing  the  two8.  The  transition  strip  is 
1120  grit  0.125"  wide  at  the  leading  edge.  In  these  figures  "dynamic 


separation"  implies  the  flow  varies  between  attached  and  separated  in  an 
intermittent  way  as  shown  in  Figure  11  (which  is  also  taken  from  Ref.  8). 
The  attached  Interval  varies  from  2-30  ellipse  revolutions,  as  is  also  true 
in  the  separated  state. 

The  data  allows  comparison  of  separated  and  attached  flow  for  the 
same  a,  k  and  Reynolds  number.  Figure  12  shows  the  first  harmonic  of 
Cp  as  a  function  of  x/c  and  Figure  13  shows  the  phase  lag  under  the 
same  conditions.  The  non- attached  unsteady  flow  seems  to  be  convected 
over  the  front  40%  of  the  chord,  and  acts  as  if  it  were  attached  over 
the  last  60%.  This  supposition  is  based  upon  a  comparison  of  and 

its  phase  in  attached  flow.  At  the  present  time  it  is  not  clear  what  is 
being  convected  at  approximately  0.55  Ua,  in  the  separated  flow 
regime.  However  the  data  (Fig.  14)  suggests  that  whatever  it  is,  it  is 
roughly  the  same  whether  the  boundary  layer  has  natural  or  artificial 
transition.  Figure  15  shows  the  variation  of  convected  velocity  with 
"k".  This  variation  is  consistent  with  that  reported  by  Carta  et 
el.*3.  we  note  too  that  when  Ap  is  added  to  AP»  t*1®  angle  of  attack 
at  stall  is  increased  by  the  unsteady  process,  as  reported  by  McCrosky 
et  al.13. 

We  note  the  oscillating  variation  of  term  Cp  on  the  upper  surface 
over  the  front  half  of  the  airfoil  acts  as  a  convected  process,  not  an 
acoustic  process,  because  a  change  in  speed  of  about  50%  does  not  change 
the  behaviour  (Fig.  16). 

An  examination  of  two  summary  plots  (Figs.  17A,  B)  for  a  range  of 
"k's",  at  a«15  degrees  and  for  natural  transition  shows  the  amplitude 
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and  phase  of  the  fundamental  of  the  unsteady  pressure  coefficient  (Cp) 
varies  from  the  leading  edge  to  the  trailing  edge  as  follows: 


k  thumps  ivalleys  phase  speed 


comments 


0.5  1  1 

1.0  2  2 

2.0  3  2 

3.9  3  2(3) 

6.4  4  4 


slow  downstream,  sign 
change  at  .8.  Last 
value  of  the  slope  is 
slightly  negative 
(-.08) 

.14  over  last  70% 
of  chord  0  over 
1st  30% 


0  over  1st  20% 

.27  over  last  80% 


.45  over  1st  50% 
zero  after 


.55  over  1st  40% 
zero  or  slightly 
negative  last  60% 


rises  slowly,  sharp 
drop  off  at  TE 


more  or  less  like 
k-.5  for  front  1/2 
chord,  oscillation 
amplitude  grows 
toward  TE 

more  or  less  const 
ampl,  rising  slope 
toward  TE 

decreasing  ampl. 
towards  TE  flat 
at  TE 

decreasing  ampl. 
toward  TE,  rising 
characteristic  at 
TE,  flat  at  midchord 


Three  features  stand  out  on  this  chart.  The  number  of  cycles 
increase  with  "k".  The  character  of  the  distribution  changes 
drastically  with  "k".  For  these  conditions,  at  low  "k"  the  amplitude 
increases  with  "k",  strongly  so,  over  the  last  20%  of  the  airfoil. 
There  is  little  going  on  over  the  front  of  the  airfoil.  At  higher  "k" 
the  oscillation  is  maximum  near  the  leading  edge  and  decreases  as  one 
moves  aft  towards  the  trailing  edge.  Further  the  convection 


characteristic  is  the  strongest  where  the  fundamental  amplitude  of  the 


unsteady  pressure  is  the  largest.  In  the  case  of  artificial  transition, 
the  indication  is  a  change-over  from  trailing  edge  separation  at  low 
"k"  to  a  leading  edge  separation  at  higher  "k".  This  dependence  at  high 
"k"  is  extraordinary  because  at  high  "k"  the  mean  plus  unsteady  pressure 
gradient  effects  are  largest  at  the  trailing  edge.  Mote  the  unsteady 
flow  field  seems  to  be  superimposed  on  the  relatively  constant  pressure 
mean  flow  field  on  the  upper  surface. 

Clearly  this  flowfield  is  not  understood. 

Comments  on  Trailing  Edge  Conditions 

The  data  from  References  1  and  2  show  the  pressure  difference 
between  the  upper  and  lower  surface  seems  to  approach  zero  as  the 
trailing  edge  is  approached.  At  low  values  of  k  the  variation  of  Ap  is 
approximately  as  the  square  root  of  distance  from  the  trailing  edge.  As 
k  increases  the  variation  is  more  nearly  linear.  Examination  of  Figures 
7  and  8  indicates  the  pressure  difference  is  finite  and  varies  smoothly 
as  the  trailing  edge  is  approached,  whether  the  flow  is  separated  or 
not.  This  implies  the  extended  Kutta  condition14 


1,2  2, 

—  (u  -  u  ) 
2  u  l 


dr 

dt 


may  well  be  applicable  here.  (Note  here  uu  is  the  tangential 


velocity,  just  outside  the  boundary  layer  at  the  upper  separation 
point,  while  u^is  the  tangential  velocity  just  outside  the  boundary 
layer  at  the  separation  point  on  the  lower  surface.  The  latter  point 


is  nominally  at  the  trailing  edge. )  For  rough  purposes  one  could  argue 


C 

P 


C  +  2 
P 


u  u  - 
l  l 


u  u 
u  u 


+ 


-2 


u 


L 


~2 

u 

u 


Thus  the  pressure  difference  has  a  part  proportional  to  the  mean,  one  to 
the  fundamental  frequency  and  as  indicated  by  the  third  term,  a  term 
proportional  to  the  second  harmonic,  as  shown  in  Figure  18,  for  the 
upper  surface  pressures.  However  the  fact  that  the  second  harmonic 
appears  only  near  the  trailing  edge  reduces  the  validity  of  this 
argument. 

On  the  Effect  of  Acoustic  Excitation 

The  effect  of  the  noise  generated  from  the  rotating  ellipse  was 
investigated.  This  was  accomplished  by  first  measuring  the  sound 
pressure  level  at  7  locations  near  the  airfoil.  This  pressure  was  then 
reproduced  using  a  speaker.  The  airfoil/speaker  configuration  is  shown 
in  Figure  19.  Many  difficulties  were  encountered  in  accurately 
reproducing  the  acoustic  wave  characteristics,  particularly  at  low  k. 
Therefore  only  reduced  frequencies  of  k-2 .0  and  6.4  were  reproduced. 

For  these  two  cases  the  sound  pressure  level  could  not  be  duplicated 
over  the  entire  airfoil  at  the  same  time  because  the  amplitude  decay  of 
the  pure  sound  wave  excitation  with  distance  was  different  than  the 
combination  of  effects  due  to  the  rotating  ellipse.  Therefore  three 
separate  sound  pressures  were  used,  one  corresponding  to  that  at  the 
leading  edge,  one  to  that  at  the  trailing  edge,  and  one  that  is 
somewhere  in  between.  (This  adds  a  large  degree  of  uncertainty  to  the 


validity  of  the  results. )  The  surface  pressure  distribution  and 
boundary  layer  profiles  were  measured  and  then  compared  to  those  for 
the  airfoil/ellipse  combination. 


The  conclusion-  vas  reached  that  at  low  values  of  "k"  (2  or  less) 
the  acoustic  part  of  the  response  is  some  30  decibels  below  the  measured 
pressure  with  wind  on.  At  k  of  6.4  the  acoustic  pressure  level 
recorded  on  the  airfoil  is  5-6  db  lass  than  the  wind-on  measurement. 
Hence  the  threshold  for  effects  of  compressibility  on  the  flow  must  lie 
in  the  neighborhood  of  these  higher  values  of  k.  For  the  most  part  the 
acoustic  excitation  had  a  negligible  effect  on  the  measured  boundary 
layer  characteristics. 

On  the  Measurement  of  Shear 

The  detailed  measurement  of  the  velocity  profile  shows  that  (in 
inner  coordinates)  there  is  a  region  where 


u  v 

T 


In  this  region  then  one  can  measure  u  near  the  wall  and  find,  if  v  and 
y  are  known 

2  vu  1  2 

u  »  “ —  »  —  0  C 
T  2  •  f 

y  * 

For  our  boundary  layers  if  we  can  measure  u  at  y  corresponds  to  y+  •  6 
or  less  we  can  deduce  Cf.  Such  probes  are  under  construction  and  will 


On  the  Relation  of  These  Results  to  Other  Results 


Certain  of  these  results  discussed  above  indicate  for  example  that 
the  unsteady  Reynolds  stress  term  (u'v1)  seems  to  be  present  in  all  the 
data^» 1 8  in  the  literature,  whether  the  data  was  obtained  in  a  channel 
with  forced  excitation,  or  on  an  airfoil  which  is  the  configuration  of 
our  experiment.  Lorber's  geometric  argument3  seems  to  apply  to  both 
channel  flow  and  our  external  flow.  If  variables  such  as  <u>  and  <u'v'> 
depend  only  on  distance  and  some  thickness  6(t),  such  that 
<u(y ,t)>  -  <u(y/6(t))> 

<u'v'(y,t)>  ■  <u'v' (y/fi (t) )>  ,  then 
time  and  space  derivatives  are  related  by 

1<2>  -  1<H>  .  iii  and 

3t  ay  4  dt 

a^'v'x  m  3<u* v1  >  _  jr  dg. 

at  ay  a  dt 


If  the  mean  and  first  harmonic  are  considered,  so  that 
<u>  (y,t)  -  u(y)  +  u(y )exp( i(ut-$u(y ) ) ) ,  and 
<u'v'>  (y,t)  '»  u'v'  (y)  +  uf’v'  (y)exp(i(ut-$uivi  (y) ) ) ,  time 
differentiation  and  division  produce  the  relation 


3  <u » v 1  > 

_ H _ 

a  <u> 
ay 


exp(i(4iu  (y)  -  4>u,v,(y)))  . 


Finally,  for  small  unsteady  amplitudes  and  for  profiles  similar  at  each 
time. 


—  <u'v'>  -  —  u'v'  and  —  <u>  »  —  u  ,  and  therefore 

ay  ay  ay  ay 


axp(iA$ 


u-u' v 


This  relation  predicts  the  amplitude  of  <u'v'>  based  upon  (i,  u  and 
u'v'.  Also,  since  the  left  hand  side  is  positive  and  real,  <u'v'>  and 
<u>  must  be  in  phase  if  au'v'/ay  has  the  same  sign  as  3u/3y  and 
180°  out  of  phase  if  the  signs  are  reversed.  Figure  20  shows  a 
comparison  of  computations  from  this  model  and  data  from  our  experiment, 
References  3  and  16  provide  further  details  and  comparisons  with  other 
data,  both  for  phase  and  amplitude. 

A  second  issue  here  is  the  relation  of  the  "Kutta  Condition"  in 
steady  and  unsteady  flow.  There  are  two  characteristics  of  the  airfoil 
associated  with  the  Kutta  condition.  The  first  is  the  statement  that 
the  pressure  jump  is  non-infinite  at  the  trailing  edge.  The  data  which 
has  been  developed  in  this  program  shows 

(1)  for  all  attached  flows  Ap*0  at  the  trailing  edge; 

(2)  for  all  separated  flow  &p  is  not  infinite  at  the  trailing  edge, 
This  aspect  of  the  "Kutta  Condition"  is  absolutely  satisfied. in  every 
case.  The  second  consequence  of  the  "Kutta  Condition"  is  to  fix  the 
value  of  the  circulations  about  the  airfoil.  These  calculations  for 


lamina  suggest  a  value  for  the  circulation  is 


This  value  increases  linearly  with  the  thickness  ratio  for  an 
airfoil  with  a  cusped  trailing  edge,  and  decreases  with  included 
trailing  edge  angle,  and  decreasing  Reynolds  number.  For  an  0012 
airfoil  the  net  result  is  that 

C  -  (0.9)C 

IDEAL 

where  C  is  the  zero  thickness  lamina  value  (2ira) 

"ideal 

Figure  21  shows  the  value  of  the  mean  CL  at  ot-10°  for  two 
Reynolds  numbers  as  a  function  of  k.  Two  clear  conclusions  are  reached. 
First  CL  mean  increases  with  Reynolds  number,  and  second  CL 
increases  with  reduced  frequency.  Figure  22  shows  the  unsteady 
increment  is  of  the  order  of  0.15.  This  represents  a  substantial 
improvement  in  effectiveness  of  developing  lift  on  heavily  loaded 
airfoils  in  unsteady  flow.  The  improvement  in  the  ability  of  an  airfoil 
to  generate  lift  in  unsteady  flow  is  in  substantial  agreement  with  the 
data  in  Ref.  13.  This  improvement  comes  about  in  spite  of  the  great 
difference  in  the  source  of  the  unsteady  flow. 

These  sorts  of  results  for  the  lift  coefficient  were  not 
unexpected.  However  the  data  on  the  unsteady  pressure  coefficient  in 
Figs.  12-17  have  yet  to  be  explained.  The  variation  of  Cp  with  x/C 
shown  there  is  unexpected  and  unexplained.  (These  matters  are  under 
study1 5.)  This  confusing  situation  is  further  complicated  by 


data  (Figs.  11,  23)  which  demonstrates  intermittent  separation  in 
unsteady  flow. 

To  further  compound  the  complexity  of  the  situation,  Figure  24 
shows  a  rough  plot  of  our  boundary  layer  data  on  a  Clauser  [a/6  -  6(G)] 
plane.  The  steady  state  "self-similar"  profiles  are  shown  as  flat 
plate,  Clauser  #1  and  Clauser  #2.  The  unsteady  data  at  small  a  is  well 
away  from  this  line  (although  the  steady  data  are  not  too  far  from  it) 
lying  at  a  value  of  "G"  that  is  too  high.  That  is,  the  unsteady  losses 
exceed  those  expected  in  steady  flow.  As  k  increases,  the  relative  size 
of  the  unsteady  losses  decrease.  At  higher  a,  the  unsteady  losses  seem 
to  be  too  small  compared  to  the  "self-similar"  thickness.  The 
implications  here  seem  simple.  As  the  Stokes  layer  becomes  thinner  with 
respect  to  the  boundary  layer  thickness  (i.e.,  as  k  increases)  more  and 
more  of  the  boundary  layer  acts  approximately  as  a  quasi-steady  layer. 
The  smaller  unsteady  amplitudes  seen  at  higher  frequencies  may  also  be  a 
factor.  The  consequences  of  this  are  treated  in  Ref.  16. 

Closing  Remarks 

At  the  present  time,  many  of  the  characteristics  of  the  unsteady 
flow  field  as  measured  in  this  program  are  unexplained.  These  vary  from 
the  frequency  doubling  at  the  trailing  edge  for  low  frequency  separated 
flows  to  the  appearance  of  a  transition  from  a  downstream  convected 
process  to  a  standing  wave  as  one  proceeds  downstream  following  a 
leading  edge  separation.  Understanding  of  some  of  these  details  may  or 


may  not  be  improved  by  measurement  of  shear  at  the  wall  under  these 
unsteady  flow  conditions.  In  other  words  the  primary  question  is  which 


of  these  characteristics  ara  a  result  of  this  particular  test 
configuration  and  which  are  characteristic  of  general  unsteady  flows. 

At  the  present  time  the  results  suggest  that  the  determination  of  the 
characteristics  of  the  separation  point  may  be  a  clarifying  result. 
Because  of  the  detailed  balance  at  the  separation  point,  it  is  suggested 
that  non-obtrusive  instrumentation  may  be  the  only  way  one  can  obtain 
unbiased  data  under  these  circumstances. 

Status  of  Research 

We  have  met  the  following  program  goals  during  the  last  grant  year. 

( 1 )  The  unsteady  boundary  layer  computer  program  is  operational  using 
an  algebraic  turbulence  modal.  Work  is  continuing  on  a  version  using  a 
three  transport  equation  turbulence  model. 

(2)  We  have  demonstrated  the  two-dimensional  character  of  the  flow  and 
find  it  satisfactory,  except  possibly  near  the  leading  edge  laminar 
separation  bubble  or  at  the  separation  point  if  transition  is  forced. 

(3)  The  behavior  of  the  unsteady  turbulent  near  wake  has  been  studied 
for  a  wide  range  of  flow  conditions,  revealing  many  interesting 
features. 

(4)  In  this  experiment  at  low  reduced  frequencies  the  unsteady  flow  is 
primarily  an  upwash  effect,  and  the  flow  behavior  is  like  that 
predicted  by  unsteady  potential  flow  of  an  incompressible  fluid.  At 
higher  reduced  frequencies  the  pressure  on  acoustic  wave  may  account 
for  up  to  10%  of  the  excitation.  The  boundary  layer  characteristics 
seem  to  be  unaffected  by  the  acoustic  excitation. 

(5)  More  of  the  test  matrix  in  a,  k  and  R,  has  been  filled  in  for 


8 


turbulent  boundary  layer  profiles,  including  measurements  of  <u>,  <v>, 
<u'u'>,  <v'v*>  and  <u'v'>. 

(6)  The  unsteady  separated  flow  zone  has  been  mapped  out.  Although  the 
surface  pressure  distribution  has  revealed  many  interesting 
characteristics,  the  details  are  not  yet  completed.  (For  example,  we 
have  yet  to  measure  the  motion  of  the  unsteady  separation  point. ) 
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RELATION  OF  OUR  DATA  TO 
CLAUSER'S  EQUILIBRIUM 


FIGURE  24s  RELATION  OF  MEAN  PROFILES  TO  SIMILARITY  SOLUTION 


